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This paper presented a multi-objective optimization of a combined cooling, heating and power system 
(CCHP) driven by solar energy. The flat-plate solar collector was employed to collect the solar radiation 
and to transform it into thermal energy. The thermal storage unit was installed to storage the thermal 
energy collected by the collectors to ensure a continuous energy supplement when solar energy was 
weak or insufficient. The CCHP system combined an organic Rankine cycle with an ejector refrigeration 
cycle to yield electricity and cold capacity to users. In order to conduct the optimization, the 
mathematical model of the solar-powered CCHP system was established. Owing to the limitation of 
the single-objective optimization, the multi-objective optimization of the system was carried out. Four 
key parameters, namely turbine inlet temperature, turbine inlet pressure, condensation temperature 
and pinch temperature difference in vapor generator, were selected as the decision variables to examine 
the performance of the overall system. Two objective functions, namely the average useful output and the 
total heat transfer area, were selected to maximize the average useful output and to minimize the total 
heat transfer area under the given conditions. NSGA-II (Non-dominated Sort Genetic Algorithm-II) was 
employed to achieve the final solutions in the multi-objective optimization of the system operating in 
three modes, namely power mode, combined heat and power (CHP) mode, and combined cooling and 
power (CCP) mode. For the power mode, the optimum average useful output and total heat transfer area 
were 6.40 kW and 46.16 m 2 . For the CCP mode, the optimum average useful output and total heat trans¬ 
fer area were 5.84 kW and 58.74 m 2 . For the CHP mode, the optimum average useful output and total heat 
transfer area were 8.89 kW and 38.78 m 2 . Results also indicated that the multi-objective optimization 
provided a more comprehensive solution set so that the optimum performance could be achieved accord¬ 
ing to different requirements for system. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Nowadays, energy consumption and carbon emission are 
becoming increasingly urgent issues. The multi-demand of energy 
for users has led to plenty efforts made by scientists to integrate 
multi-energy systems to achieve a higher efficiency. In this situa¬ 
tion, the combined cooling, heating and power system (CCHP) 
has proved to be an efficient way to achieve energy saving and eco¬ 
nomic saving as well as to reduce the emission of pollutants [1,2], 
Absorption chiller, compression refrigerator or ejector chiller can 
be used in the CCHP system to produce cooling capacity. 

A large amount of works focused on CCHP systems have been 
carried out on the performance optimization considering the cost 
criterion. The linear programming model and the non-line 
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programming model were employed to conduct the economic 
optimization of the CCHP system to obtain the minimum value of 
the cost of energy [3-5]. Ren et al. [6] developed a mixed integer 
nonlinear programming model to optimize the annual cost of a res¬ 
idential CCHP system. Tichi et al. [7] employed a particle swarm 
algorithm to minimize the cost function for different CHP and 
CCHP systems in an industrial dairy unit. It could be confirmed that 
the CCHP system was more efficient in energy saving. 

Except for the economic index, the CCHP system performance 
can be additionally determined by thermodynamic and environ¬ 
mental criteria. Several articles also carried out additional analysis 
from the perspective of thermodynamics or environmental aspects. 
Chen et al. [8] conducted an exergoeconomic performance optimi¬ 
zation of a CCHP system with a closed Brayton cycle. The fossil fuel 
energy savings and the C0 2 emissions were compared with the 
optimization results of a conventional system to evaluate the sys¬ 
tem performance more comprehensively. Oh and Seo et al. [9,10] 

















290 


M. Wang et al./Energy Conversion and Management 89 (2015) 289-297 


Nomenclature 



A 

area, m 2 

Subscripts 

Bo 

boiling number 

a 

ambient; absorption 

c P 

specific heat at constant pressure, kj/kg K 

b 

beam radiation 

D 

diameter, m 

c 

cold side 

E 

exergy, kW 

cond 

condenser 

F 

collector efficiency factor 

cool 

cooling 

Fr 

heat removal factor 

CND 

condenser 

h 

convection heat transfer coefficient, W/m 2 K; enthalpy, 

d 

diffuse radiation; daily 


kj/kg 

Eq 

equivalent 

I 

hourly radiation, W/m 2 

EVP 

evaporator 

m 

mass flow rate, kg/s 

h 

hot side; hydraulic diameter of flow channel 

Nu 

Nusselt number 

heat 

heat 

P 

pressure, MPa 

heater 

heater 

Pr 

Prandtl number 

in 

input 

Q 

heat rate, kW 

i 

inlet; inner 

R 

tilt factor for radiation 

i 

loss; liquid phase 

Re 

Reynolds number 

L 

liquid in the water tank 

S 

incident solar flux, W/m 2 

FC 

flat-plate collector 

T 

temperature, I<; Time, h 

load 

load 

U 

loss coefficient, W/m 2 K; heat transfer coefficient, 

m 

mean 


W/m 2 K 

max 

maximum 

V 

volume, m 3 

min 

minimum 

W 

pitch of tube, m; power, kW 

net 

net 

X 

vapor quality 

out 

output 



0 

outer 



p 

collector plate 

Greek letters 

pinch 

pinch point 

a 

absorptivity 

pump 

r 

pump 

reflected radiation 

a 

Chevron angle of the plates 

s 

isentropic point; Sun 

y 

ratio of the average useful output to the total heat trans- 

solar 

solar 


fer area 

t 

tank 

6 

thickness, cm 

tot 

total 

Atm 

log mean temperature difference between hot side and 

turbine 

turbine inlet 


cold side, K 

TBN 

turbine 

'/ 

efficiency, % 

u 

useful heat gain 

A 

P 

thermal conductivity, W/m I< 
density, kg/m 3 

V 

VG 

vapor phase 
vapor generator 

T 

transmissivity; Time, h 

w 

water in the tank 

<t> 

latitude, ° 

waterl 

water into the solar collector 



water2 

water into the vapor generator 


used the mixed integer linear programming models to minimize 
the total annual cost of trigeneration systems in commercial and 
residential sectors. Arcur et al. [11] conducted the economic opti¬ 
mization of a trigeneration system in an energetically complex 
context by using a mixed linear programming model. Then, the 
energetics and environmental analysis were conducted based on 
results of the economic optimization to achieve a better system 
performance. 

In order to evaluate the system performance comprehensively, 
certain researches focused on the optimization of CCHP systems 
were performed concerning simultaneously energetic, economic 
and environmental aspects. The single-objective optimizations 
for the CCHP system [12-15] were conducted by integrating pri¬ 
mary energy consumption (PEC), operation cost (OC), and carbon 
dioxide emissions (CDE) or other index as the objective function 
respectively. The results demonstrated that there was not an opti¬ 
mal condition where each objective all reached the optimum value. 
Liu et al. 16] utilized the sequential quadratic programming (SQP) 
algorithm to solve the optimization problem of a CCHP system. The 
ratio of electric cooling to cool load was varied to minimize the 


objective function which was made up of three parts, namely the 
primary energy consumption, the hourly and annually total cost 
and the greenhouse gas (GHG) emission. 

However, much work has been done on the single-objective 
optimization of the CCHP system, while few studies have been 
investigated in the multi-objective optimization. Abdollahi and 
Meratizaman [17-19] used the genetic algorithm to conduct the 
multi-objective optimization of a CCHP system. The objectives 
functions were selected from the view of thermodynamic, eco¬ 
nomic and environmental aspects. A two-stage optimal planning 
and design method was applied for a CCHP system [20]. Non-dom- 
inated Sorting Genetic Algorithm-II (NSGA-II) was employed in 
multi-objective optimization on the first stage, and the mixed-inte¬ 
ger linear programming algorithm on the second stage. 

All of these researches above have been devoted to CCHP pow¬ 
ered by fossil fuels which will be limited with the growing crisis of 
energy shortage and are harm to environment with CO x and NO x 
emissions. The renewable energies, including hydraulic energy, 
hydrogen energy, solar energy and so on, were clean and 
inexhaustible, which would contribute to resolve the problem of 
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energy crisis and emission of pollutions. In this situation, integrat¬ 
ing the renewable energy technology to the CCHP system seems to 
be an efficient way to improve the energy efficiency and to reduce 
the pollutant emission to the environment. Solar energy, as a kind 
of renewable energy resources, is widely applied to the trigenera¬ 
tion system [21 . Compared with the independent cooling, heating 
and power systems powered by solar energy, the solar-powered 
CCHP system has attracted more attention of researchers. Because 
in the combined configuration, the three subsystems can share 
same solar collection and storage system, pipes or auxiliary 
equipment, while the independent systems require respective 
independent equipment which may cause a higher cost. Wang 
et al. [22] combined a Brayton cycle and a transcritical C0 2 refrig¬ 
eration cycle to perform the CCHP system driven by solar energy. 
Furthermore, they [231 also examined a solar-powered CCHP sys¬ 
tem for a specific building based on organic working fluid. Para¬ 
metric analysis was conducted to track the operation 
performance and to analyze the influence of key parameters on 
the performance from both thermodynamic and economic aspects. 

From the view of engineering application, it is necessary to con¬ 
duct the optimization of solar-powered CCHP which is helpful to 
the design of system. Wang et al. [24,25] proposed a new CCHP 
system powered by solar energy. Genetic algorithm was employed 
to carry out the system optimization with modified system effi¬ 
ciency or the exergy efficiency as the objective function respec¬ 
tively. Tora and El-Halwagi [26] used a nonlinear programming 
formulation to minimize the total annualized cost of a trigenera¬ 
tion system driven by solar energy and fossil fuels. Since there 
are many factors having the significant effects on the system per¬ 
formance from thermodynamic, economic and environmental 
aspects, the single-objective optimization could not evaluate the 
performance of the solar-powered CCHP system comprehensively. 
Multi-objective optimization should be urgent to be conducted. 

In this study, the multi-objective optimization of a combined 
cooling, heating and power system (CCHP) driven by solar energy 
is conduct to obtain the optimum performance from the view of 
two aspects. Non-dominated Sorting Genetic Algorithm-11 is 
employed to conduct the multi-objective optimization. The trigen¬ 
eration system is self-sufficient by using flat-plate collectors, meet¬ 
ing heating, cooling and power demands of users simultaneously. A 
thermal storage tank is employed to store solar energy collected by 
collectors and to supply energy to the CCHP subsystem continu¬ 
ously even when the solar radiation is insufficient. The average 
useful output and the total heat transfer area are selected as the 


objective functions. Turbine inlet temperature, turbine inlet pres¬ 
sure, condensation temperature and pinch temperature difference 
in vapor generator are selected as the decision variables. 

2. System description 

The CCHP system powered by solar energy consists of two sub¬ 
systems, namely the solar collection subsystem and the CCHP sub¬ 
system, as shown in Fig. 1. 

The solar collection subsystem is made up of a solar collector, a 
thermal storage tank and an auxiliary heater. The flat-plate collec¬ 
tor is selected to collect solar radiation due to its low cost and wide 
application. A thermal storage tank is employed to correct the mis¬ 
match between the supply of the solar energy and the demand of 
thermal source consumed by the CCHP subsystem, thus the system 
could operate continuously and stably. On rainy days or at night, 
when the solar radiation is insufficient and the temperature of 
thermal storage tank drops below the allowable reference temper¬ 
ature, the auxiliary heater is activated to elevate the outlet temper¬ 
ature of thermal storage tank. The heat-transfer medium in the 
solar collection subsystem is water for its low cost and large heat 
capacity. 

The CCHP subsystem, which mainly combines an organic 
Rankine cycle (ORC) with an ejector refrigeration cycle, consists 
of a vapor generator, a turbine, an evaporator, a heater, a 
condenser, a recuperator, a throttle valve, an ejector, a pump, 
and several regulation valves. In the vapor generator, the liquid 
working fluid with high pressure is vaporized by absorbing heat 
from the hot water, and then the high pressure vapor is divided 
into three streams. One stream is delivered to a turbine where it 
expands to a low pressure to drive a generator to produce electric¬ 
ity. Another stream flows into an ejector as a primary flow to suck 
the secondary flow from an evaporator with very low pressure. 
Afterwards, the primary flow and the secondary flow are mixed 
together to a mid-pressure leaving the ejector. The third stream 
is delivered to a heater to provide heat energy to users and then 
flows into a recuperator. The stream leaving the ejector is mixed 
with the turbine exhaust and then condensed to liquid state in a 
condenser by rejecting heat to the environment. The liquid work¬ 
ing fluid is also divided into two streams. One stream is throttling 
down to a very low pressure across a throttle valve and then is 
vaporized in an evaporator by absorbing heat to produce cooling 
effect. The other stream is pumped to a high pressure and enters 
a recuperator where it is mixed with the stream out of the heater. 



Fig. 1 . Schematic diagram of the solar-powered CCHP system. 
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Finally, the working fluid leaving the recuperator is delivered back 
to the vapor generator. 

In this case, the CCHP system can operate in three modes, 
namely power mode, combined heat and power (CHP) mode, and 
combined cooling and power (CCP) mode. During a year-round, 
power load is needed by the users all the time. In spring or autumn, 
heat and cooling is not necessary so that the system operates in 
power mode by closing the regulation valve 1 and 2. In winter, reg¬ 
ulation valve 1 is opened to make the system operate in CHP mode. 
In summer, regulation valve 2 is opened while the regulation valve 
1 is closed so that the system could discharge cooling and power to 
load in CCP mode. 

As to the working fluid in the CCHP subsystem, we consider 
thermodynamic properties and environmental friendliness, and 
finally select R245fa as the working fluids of the system due to 
its good performance and little negative impact on environment. 

3. Mathematical models and performance criteria 

3.1. Mathematical models 


The mathematical models of the solar-powered CCHP system 
mainly consist of the solar collection models and the CCHP models. 

The flat-plate collector is selected to absorb the solar radiation. 
The hourly radiation falling on a tilted surface can be given by 

he = IbRb + + (fb + h)F r (1) 


where I b and J d are the hourly beam and diffuse radiation that the 
collector receives, R bl R d , and R,- are defined as tilt factors for differ¬ 
ent radiations. The solar radiation and the ambient temperature are 
assumed to be constant within one hour. 

The incident solar flux absorbed by the absorber plate is calcu¬ 
lated by 

S = I b Kb(ra) b + [J d R d + (J b + / d )R r ](Toc) d (2) 


where (ta) b and (ra) d represent the transmissivity-absorptivity 
product falling on the collector. 

Considering the heat loss in the collector, the useful heat gain 
rate of the collectors can be achieved by introducing the collector 
heat-removal factor F R [27] 


F r = 


ftlwaterl 6p 
UtA p 


, , F'U t A p 


Q U = F R A P .[S-U,(T L -T a ) 


(3) 

(4) 


where the collector efficiency F is determined by [28] 


WU, 1_1_u Js_ j_L_ l 

1 l^KW-DW+Do] ^ fc a D 0 ^ JlDj/i, J 

A water storage tank is installed to act as a buffer between the 
solar collector and the bottom CCHP subsystem. The tank is 
assumed to be insulated and the water is assumed to be well- 
mixed so that the water temperature T L varies only with time. 
The following equation can be obtained from the energy balance 
in the tank [28] 

[(pVC v )„ + (pVCp)J ^f=Q.u- Qload - m\ ■ (Tl - Fa) (6) 

Furthermore, Qioad. which represents the energy discharged to 
CCHP subsystem, can be calculated as 

Qioad = m water2 Cp(T L - Ti) (7) 

The bottom CCHP subsystem combines an organic Rankine 
cycle with an ejector refrigeration cycle, yielding heat, power and 
cooling to users. To simplify the simulation of the system, several 
assumptions about the CCHP subsystem are made [24]: 


(1) The pressure drop in pipes or other components, the heat 
loss to the environment in the heat exchangers, and the fric¬ 
tional and mixing losses in the ejector are neglected. 

(2) The working fluid at the condenser outlet is saturated liquid. 
The condensation temperature is approximately 5 °C higher 
than the environment temperature. 

(3) The flow through the throttle valve is isenthalpic and the 
flow in ejector is assumed to be one-dimensional. 

For the vapor generator, the heat balance equation is: 

C p m wa ter 2 (F L - Ti) = m 3 (h 3 - h 2 ) (8) 

For the turbine and the pump, the isentropic efficiency is con¬ 
sidered. The power generated by the turbine and consumed by 
the pump are the assumed to be the power transferred from the 
fluid to the moving parts of turbine and to the fluid from the mov¬ 
ing parts of pump respectively. 

h 3 -h 4 

Wtbn = tn 4 (h 3 - h 4 ) 

_ h t , - h 6 

'/pump fti _ he 

Wpump = m,(h 1 -he) 

Therefore, the net power output can be given by 

Wnet = Wtbn - W pllmp (13) 

For the heater and evaporator, the heat and cooling output to 
users respectively are: 

Qheat = m,o • (h 3 — hjo) (14) 

Qcool = tn s ■ (hs - h 7 ) (15) 

All the heat exchangers in this system are plate heat exchanger 
(PHE) type for its high efficiency and compact structure. The total 
heat transfer area of the CCHP system is determined from a sum¬ 
mation of individual components: 

Aot = -AcnD + Ay G + A EV p + Aneater (16) 

In power mode, the heater and the evaporator do not work, so 
the heat transfer area of the heater and the evaporator are not con¬ 
sidered in the calculation of the total heat transfer area. While in 
CCP mode the heat transfer area of the heater is not considered, 
and in CHP mode the heat transfer area of the evaporator is not 
considered. 

In the heat exchanger, the working fluid may work in different 
thermodynamic states, namely superheated state, two-phase state 
and subcooled state, resulting in that the heat transfer process in 
the heat exchanger can be divided into different regions. In the 
condenser, the working fluid is liquefied from superheated state 
to saturated liquid state undergoing the superheated region and 
the two-phase region. In the vapor generator and the heater, the 
working fluid respectively works in subcooled region, two-phase 
region and superheated region. In the vapor generator, the pinch 
temperature difference is an important parameter, which is the 
minimum temperature difference between the heat source and 
the saturation temperature of working fluid at the evaporation 
pressure. In the evaporator, however, the working fluid is vapor¬ 
ized only in two-phase region. Compared to the two-phase region, 
the superheated region and the subcooled region are single-phase 
region. 

In the heat exchanger, the general formula of the heat rate in 
any region can be given by [29] 

Q. = UAAt m (17) 

The log mean temperature difference is expressed as 


(9) 

( 10 ) 

(ID 

( 12 ) 
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* . _ Afmax Afmin 

m_ ln(At max /At min ) 

The overall heat transfer coefficient is given by 
11 5 1 

u ~ hu + Am + h c 


(18) 


(19) 


where the convection heat transfer coefficient for single-phase flow 
is expressed as 


h = 


X ■ Nu 

~dT 


( 20 ) 


For the single phase region, the Nusselt number is calculated by 
using the Chisholm and Wanniarachchi correlation 


/ co\ 0.646 

Nu = 0.724 m Re°' 583 Pr 1/3 


( 21 ) 


where /I stands for the chevron angle of the plates, being expressed 
by Ref. [30], 

For the two-phase region, the heat transfer process is divided 
into relatively small sections so that the properties in each section 
can be assumed to be constant with such slight property variations. 
The Nusselt numbers are different when the working fluid is vapor¬ 
ized or condensed. For each section, the calculations of Nusselt 
number [31,32] in the hot side and in the cold side when the work¬ 
ing fluid is in two-phase region can be respectively given by 


Nu ih = 4.118Re° q 4 jPr ( 1/3 
Nui, c = 1.926Re° q 5 i Pr, 1/3 Bo° q 3 i 


1 Xm,i “F ^in.i 



( 22 ) 

(23) 


In addition, the one-dimensional model of the ejector is used to 
simplify the simulation process. The detailed description about the 
ejector model is given in Ref. [33], 


3.2. Performance criteria 


Since the quality of cold energy, heat energy and power energy 
are different, the sum of cold, heat and power could not reflect the 
output capacity of the whole system. Thus, two factors, namely a 
and /?, are introduced to modify the different price between cold, 
heat and power energy of the output of the system, given by 

Wout = W net + aQheat + Wcold (24) 


where a and ft vary with the price fluctuation of cooling, heating 
and power. In this paper, a = 0.5 and ft = 0.8 [34], 

The working condition of the system varies with time over a 
day. A daily average efficiency of the CCHP system is modified 
compared to the instantaneous efficiency to assess the whole-day 
performance of the system, shown as 


/Woutdr 

~ f Q u dr 


(25) 


The daily exergy efficiency is examined from the view of the 
second law of thermodynamics, shown as 


_ / W out di 

neX ~ /Bind! 


(26) 


where £ in represents the input exergy of the system. 

In this CCHP system, £ in can be calculated as 

Bin = B so |ar (27) 


where £ solar is the exergy input of the solar collector, defined as 


Bsolar — Ap/pc 


1 + 


1 

3 




(28) 


T s is the sun temperature and equals to 6000 K according to Ref. 

[35], 

In order to assess the system performance from both thermody¬ 
namic and economic aspects, a thermo-economic index y, named 
as output-to-area ratio, is defined as the ratio of the average useful 
output to the total heat transfer area of the CCHP system in the sin¬ 
gle-objective optimization, calculated by 


where W out m presents the average useful output of the CCHP sys¬ 
tem over a day 


kPout.m — 


\! 


W out dT 


(30) 


4. Results and discussions 

In order to achieve the optimal performance of the solar-pow¬ 
ered CCHP system, we conduct the optimization in the three typi¬ 
cal conditions year-round. The Genetic Algorithm (GA) is employed 
to find the optimum performance of system. Four key thermody¬ 
namic parameters, namely turbine inlet temperature, turbine inlet 
pressure and condensation temperature and pinch temperature 
difference are chosen as the decision variables. The main thermo¬ 
dynamic parameters for the simulation of the solar-powered CCHP 
system are listed in Table 1. 

The ranges of the decision variables in the optimization are 
listed in Table 2. The thermodynamic properties of the working flu¬ 
ids are calculated by REFPROP 9.01 [36] developed by the National 
Institute of Standards and Technology of the United States. 

4.1. Single-objective optimization 

For the single-objective optimization, the objective function is 
the defined thermo-economic index, y. The optimum system per¬ 
formance and the corresponding combination of the decision vari¬ 
ables are shown in Table 3. The results demonstrate that the 
optimum turbine inlet temperature in power, CCP and CHP mode 
are 77.59 °C, 78.01 °C and 67.63 °C respectively, and the turbine 
inlet pressure is approximately equal to the corresponding satu¬ 
rated pressure. Because the cooling and power are both high qual¬ 
ity energies, generating power or cooling energy needs to consume 
more energy. Therefore, the optimal turbine inlet temperatures of 


Table 1 

Simulation conditions for the solar-powered CCHP system. 


Term 

Value 

Simulation date in winter 

December 21st 

Simulation date in summer 

July 20th 

Simulation date in spring/autumn 

March 21st 

Collector tilt angle (CHP mode)/ 0 

61 

Collector tilt angle (CCP mode)/ 0 

7 

Collector tilt angle (power mode)/ 0 

35 

Outer diameter of the absorber tube/m 

0.018 

Inner diameter of the absorber tube/m 

0.014 

Surface area of one absorber plate/m 2 

1.76 

Number of absorber plates 

500 

Number of water tank 

20 

Evaporator temperature/°C 

5 

Heater outlet temperature/°C 

50 

Turbine isentropic efficiency/% 

60 

Pump isentropic efficiency/% 

70 
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Table 2 

Data of the parameter optimization. 



651 

Term 


Value 

60- 

Population size 


60 


Crossover probability 


0.8 

55- 

Mutation probability 


0.02 


Stop generation 


100 

50- 

Range of turbine inlet temperature/°C 


65-80 

£ 

Range of turbine inlet pressure/MPa 


0.40-0.80 


Range of condensation temperature/°C 

CCP mode 

35-45 

^ 45- 


Power mode 

15-25 



CHP mode 

10-20 

40- 

Range of pinch temperature difference/°C 


5-10 



35- 


■ Power-Spring/Autumn 



Table 3 

Single-objective optimization results of the solar-powered CCHP system in three 
operation modes. 


Term 

Power mode 

CCP mode 

CHP mode 

Output-to-area ratio/kW m -2 

0.191 

0.125 

0.270 

Average useful output/kW 

6.96 

7.16 

8.45 

Total heat transfer area/m 2 

36.43 

62.76 

31.34 

Daily average efficiency/% 

9.75 

8.76 

13.81 

Daily exergy efficiency/% 

6.16 

6.16 

8.41 

Turbine inlet temperature/°C 

77.59 

78.01 

67.63 

Turbine inlet pressure/MPa 

0.743 

0.751 

0.572 

Condensation temperature/°C 

20.00 

35.00 

10.00 

Pinch temperature difference/°C 

10.00 

10.00 

9.97 


system in power mode and in CCP mode are higher than that of 
system in CHP mode. In addition, the optimum condensation tem¬ 
peratures are close to their minimum values while the optimum 
values of pinch temperature difference is close to their maximum 
values in their design ranges. The optimal system performances 
under this condition of optimization are 0.191 kW/m 2 for power 
mode, 0. 125kW/m 2 for CCP mode and 0.270 kW/m 2 for CHP 
mode, meanwhile the exergy efficiencies in power, CCP and CHP 
mode reach 6.16%, 6.16% and 8.41% respectively. 


30- 


6.0 


6.1 


6.2 


6.3 

W /kW 


6.4 


6.5 


6.6 


Fig. 2. Results of multi-objective optimization in power mode. 



4.2. Multi-objective optimization 


Fig. 3. Results of multi-objective optimization in CHP mode. 


Non-dominated Sorting Genetic Algorithm-II (NSGA-II), which 
is an effective multi-objective optimization method, is employed 
to find the solution of the multi-objective optimization. Owing to 
the fluctuation of solar energy over a day, the useful output of sys¬ 
tem varies correspondingly. Therefore, the average useful output is 
defined to assess the capacity of the system. Expect for thermody¬ 
namic performance, the economic index should be another impor¬ 
tant criterion to examine the performance of the CCHP system. 
Actually, the size of turbine is independent on the power capacity 
to some degree because larger capacity turbine can use high rota¬ 
tional speed design to decrease its size. Thus, the size of turbine is 
not sensitive to the capacity which is affected by thermodynamic 
parameters. So is the pump. But the areas of heat exchangers are 
very sensitive to thermodynamic parameters. As the heat transfer 
area is the mainly effect factor of the capital cost of the heat 
exchanger, the total heat transfer area, which consists of the heat 
transfer area of vapor generator, evaporator, condenser and heater, 
is selected as another performance criterion. Eventually, the objec¬ 
tives of the multi-objective optimization are to maximize the aver¬ 
age useful output Wout, m and to minimize the total heat transfer 
area A tot . 

Figs. 2-4 show the Pareto frontier solutions for the solar 
-powered CCHP system in three operation modes of multi-objec¬ 
tive optimization, which reveal the conflicts between average use¬ 
ful output and total heat transfer area. Comparison among the 
three figures demonstrates that the total heat transfer area in 



W a JkW 


Fig. 4. Results of multi-objective optimization in CCP mode. 


CHP mode is smaller than those in other two modes. The reason 
is that the cooling and power are both high quality energy which 
requires a larger total heat transfer area in a whole system for 
same capacity output. However, the larger average useful output 
happens to the system in CCP mode. Because the solar radiation 
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in summer is much stronger than those in other seasons, implying 
that the system receives more heat input. Results also indicate that 
in order to achieve a large average useful output, the total heat 
transfer area should be large enough as well, while the total heat 
transfer area increases rapidly when the average useful output is 
greater than the specific value. Therefore, only seeking for a large 
average useful output is not the best choice for the optimal system 
design. 

In these figures, each point located on the Pareto frontier is 
potentially an optimum solution. A larger average useful output 
requires a larger heat transfer area. The maximum average useful 
output is defined as design point A, where the total heat transfer 
area is at its maximum value as well. The minimum total heat trans¬ 
fer area is defined as design point B, with corresponding minimum 
average useful output. Design point A and B respectively imply the 
optimum solutions of single-objective function where the average 
useful output or the total heat transfer area is selected as the objec¬ 
tive function. Since it is impossible to make each objective at its 
optimum value simultaneously, the absolutely optimum solution 
of multi-objective does not exist. The point satisfying the result with 
the minimum total heat transfer and the maximum average useful 
output simultaneously does not exist in the optimal solutions, 
shown as Point D in Figs. 2-4. Therefore, in this paper, the point 
which is the closest point to the specific point D, is decided to be 
the final solution in the multi-objective optimization. Thus, to 
achieve the final solution, a process of decision-making is required 
with the aid of a hypothetical point D [37] which is the intersection 
of maximum output and maximum heat transfer area and is not 
located on the Pareto frontier. In this study, the closest point C of 
the Pareto frontier to the hypothetical point D is considered as the 
final solution. At point C, the system achieves the best possible val¬ 
ues considering two objective functions. Tables 4-9 list the opti¬ 
mum objective function values and the corresponding decision 
variables for points A-C on the Pareto frontier and the correspond¬ 
ing detailed parameter distribution on the system in power, CHP 
and CCP modes respectively. Results indicate that under the condi¬ 
tion of the optimum solutions of the multi-objective optimization of 
the combined system in CCP and CHP mode, the cooling-to-power 
ratio and heat-to-power ratio are 1.29 and 1.01 respectively. 

Fig. 5 depicts the distribution of turbine inlet pressure for the 
Pareto frontier in its design range in different modes. It can be 
observed that the optimum solutions of turbine inlet pressure in 
CCP mode and CHP mode respectively locate in the range of 
0.5-0.75 MPa, while in power mode the solutions are around 
0.7 MPa. 

Fig. 6 depicts the distribution of turbine inlet temperature for 
the Pareto frontier in its design range in different modes. It is 
shown that the optimum solutions of turbine inlet temperature 
locate in the range of 70-78 °C in CCP mode and CHP mode and 
75-77 °C in power mode respectively. The result of the optimum 
turbine inlet temperature is corresponded to that of the optimum 
turbine inlet pressure. In summer when the solar radiation is 


Table 4 

Optimum design objective values and corresponding decision variables for points A-C 
on the Pareto optimal front in power mode. 


Term 

A 

B 

C 

Average useful output/kW 

6.52 

6.15 

6.40 

Total heat transfer area/m 2 

50.60 

38.45 

46.16 

Daily average efficiency/% 

9.61 

9.61 

9.54 

Daily exergy efficiency/% 

6.52 

6.15 

6.40 

Turbine inlet temperature/°C 

75.04 

76.44 

75.79 

Turbine inlet pressure/MPa 

0.695 

0.721 

0.689 

Condensation temperature/°C 

20.00 

20.03 

20.04 

Pinch temperature difference/°C 

5.43 

9.95 

6.76 


Table 5 

Detailed parameter distribution on the system in power mode. 


State 

t/°C 

p/MPa 

p/kg m 3 

hi kj kg-' 

s/kj kg- 1 K-' 

ml kgs 1 

1 

20.37 

0.721 

1352.81 

226.53 

1.092 

0.426 

2 

20.37 

0.721 

1352.81 

226.53 

1.092 

0.426 

3 

76.41 

0.721 

40.08 

459.39 

1.778 

0.426 

4 

35.55 

0.123 

6.70 

433.65 

1.799 

0.426 

5 

35.55 

0.123 

6.70 

433.65 

1.799 

0.426 

6 

20.03 

0.123 

1351.94 

225.90 

1.091 

0.426 


Table 6 

Optimum design objective values and corresponding decision variables for points A-C 
on the Pareto optimal front in CCP mode. 


Term 

A 

B 

C 

Average useful output/kW 

7.40 

2.13 

5.84 

Total heat transfer area/m 2 

74.66 

42.07 

58.74 

Daily average efficiency/% 

8.71 

2.26 

6.90 

Daily exergy efficiency/% 

6.36 

1.83 

5.02 

Turbine inlet temperature/°C 

77.68 

70.65 

77.17 

Turbine inlet pressure/MPa 

0.744 

0.499 

0.694 

Condensation temperature/°C 

35.11 

42.36 

37.72 

Pinch temperature difference/°C 

6.90 

9.76 

9.61 


Table 7 

Detailed parameter distribution on the system in CCP mode. 

State 

t/°C 

p/MPa 

p/kg m 3 

hi kj kg-' 

s/kj kg-' K-' 

ml kgs ' 

1 

38.01 

0.694 

1304.12 

249.98 

1.170 

0.527 

2 

38.01 

0.694 

1304.12 

249.98 

1.170 

0.527 

3 

77.17 

0.694 

38.06 

460.97 

1.784 

0.527 

4 

50.52 

0.232 

12.41 

444.82 

1.797 

0.293 

5 

55.36 

0.232 

12.16 

449.62 

1.811 

0.557 

6 

37.72 

0.232 

1303.24 

249.47 

1.169 

0.557 

7 

5.00 

0.066 

18.49 

249.47 

1.178 

0.030 

8 

5.00 

0.066 

3.99 

408.12 

1.748 

0.030 

9 

60.72 

0.232 

11.91 

454.93 

1.827 

0.264 


Table S 

Optimum design objective values and corresponding decision variables for points A-C 
on the Pareto optimal front in CHP mode. 


Decision variables 

A 

B 

C 

Average useful output/kW 

9.28 

7.67 

8.89 

Total heat transfer area/m 2 

53.17 

28.90 

38.78 

Daily average efficiency/% 

13.63 

13.91 

13.59 

Daily exergy efficiency/% 

9.24 

7.63 

8.85 

Turbine inlet temperature/°C 

71.91 

76.69 

72.55 

Turbine inlet pressure/MPa 

0.802 

0.716 

0.513 

Condensation temperature/°C 

10.05 

10.26 

10.11 

Pinch temperature difference/°C 

5.14 

9.99 

8.03 


Table 9 

Detailed parameter distribution on the system in CHP mode. 


State 

t/°C 

P/MPa 

p/kg m 3 

h/kjkg ' 

s/kj kg 1 K 1 

ml kgs 1 

1 

10.34 

0.513 

1303.24 

213.42 

1.047 

0.322 

2 

15.66 

0.513 

1364.68 

220.32 

1.071 

0.371 

3 

72.55 

0.513 

27.16 

460.42 

1.799 

0.371 

4 

33.67 

0.082 

4.46 

433.21 

1.821 

0.322 

5 

33.67 

0.082 

4.46 

433.21 

1.821 

0.322 

6 

10.11 

0.082 

1378.06 

212.98 

1.047 

0.322 

10 

50.00 

0.513 

1268.18 

266.32 

1.222 

0.048 


strong, the working fluid can be heated to a higher evaporation 
temperature. 

Fig. 7 depicts the distribution of condensation temperature for 
the Pareto frontier in its design range in different modes. It is dem- 
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Fig. 5. Distribution of turbine inlet pressure for the Pareto frontier in its design 
range in different modes. 



Fig. 6. Distribution of turbine inlet temperature for the Pareto frontier in its design 
range in different modes. 



Fig. 7. Distribution of condensation temperature for the Pareto frontier in its design 
range in different modes. 


onstrated that the optimal solutions of condensation temperature 
in power mode and CHP mainly locates in around the minimum 
value of the design ranges, and the optimum solutions of turbine 
inlet pressure in CCP mode exist in the range of 35-43 °C. It can 



Fig. 8. Distribution of pinch temperature difference for the Pareto frontier in its 
design range in different modes. 

be concluded that the optimum condensation temperature is not 
always the smaller ones. The condensation temperature affects 
the heat transfer of not only condenser but also other heat 
exchangers in the system. The optimum solution would be various 
in different conditions. 

Fig. 8 depicts the distribution of pinch temperature difference 
for the Pareto frontier in different modes. The optimum pinch tem¬ 
perature difference in power mode and CHP mode locate in the 
range across the design range, meanwhile the optimum pinch tem¬ 
perature difference in CCP locates in the range of 7-10 °C. 

In general, the heat source temperature has a significant impact 
on the turbine inlet temperature, while system in CHP mode 
requires for a lower heat source temperature. Hence, most of the 
optimal turbine inlet temperature solutions for system in CHP 
mode are lower than other modes. As to the condensation temper¬ 
ature, it is obvious that the lower condensation temperature, the 
more power output, while the heat and cooling are barely affected. 
However the heat transfer area of condenser significantly increases 
when the condensation temperature is set to a low value. Finally, 
the optimal solutions of condensation temperature in power mode 
locate in the specific range. 

The turbine inlet temperature should not be designed to be too 
high or too low; meanwhile the turbine inlet pressure should 
ensure that the working fluid before entering the turbine is satu¬ 
rated or superheated vapor. The condensation temperature should 
be as low as possible in a certain range that the heat transfer area is 
in proper size. 

5. Conclusions 

In this paper, the multi-objective optimization of a solar-pow¬ 
ered CCHP system based on organic working fluid is conducted. 
The system operates in three operation modes, namely power 
mode, CHP mode and CCP mode. Turbine inlet pressure, turbine 
inlet temperature, condensation temperature and pinch tempera¬ 
ture difference are selected as the decision variables, while the 
average useful output and the total heat transfer area are selected 
as the objective functions for the optimization. NSGA-II is employed 
to achieve the optimum performance in the multi-objective 
optimization. 

The solutions of the optimization are obtained. For the power 
mode, the optimum average useful output and total heat transfer 
area are 6.40 kW and 46.16 m 2 . For the CCP mode, the optimum 
average useful output and total heat transfer area are 5.84 kW 
and 58.74 m 2 . For the CHP mode, the optimum average useful 
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output and total heat transfer area are 8.89 kW and 38.78 m 2 . 
Eventually, in the situation of this case, the turbine inlet tempera¬ 
ture should be set to be a proper value based on the requirement 
that the turbine inlet pressure could ensure that the working fluid 
in the turbine inlet is saturated or superheated vapor. On the con¬ 
dition that the heat transfer area of the condenser is not too large, 
the condensation temperature could be as low as possible. The 
pinch temperature difference can be decided in different range 
for different modes. The optimal combination of the parameters 
in single objective optimization is different from that in multi¬ 
objective optimization under the given condition. Every point of 
the solutions in the multi-objective optimization is a candidate 
of the optimum result, implying that the multi-objective optimiza¬ 
tion provides a more comprehensive solution set so that the best 
performance can be achieved according to different requirements 
for system. 
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